Irradiated Graphene Nanostructures as Nanoscopic Pressure Cells 
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High-dose irradiation of nanostructures consisting of multiple graphene shells, such as spherical 
'carbon onions' (CO) or cylindrical multi- walled carbon nanotubes (MWNT), induces shell shrink- 
age by bond reconstruction around irradiation-induced defects. This leads to build-up of internal 
stresses and to extreme pressures acting on encapsulated materials or on the carbon cores of the 
nanostructures. We formulate a model which relates the build-up of internal stresses to the point 
defect dynamics in graphene shells. Calculations are performed for the special case of irradiated 
CO. The results are in good agreement with experimental findings. 



Carbon onions (CO) consisting of concentric spher- 
ical graphene shells form under intense irradiation of 
graphitic or amorphous carbon with electrons [lH3[ or 
ions [4(. At temperatures above 600 K, in-situ anneal- 
ing of radiation defects leads to coherent graphene shells 
which contract during irradiation as manifested by de- 
creasing shell spacing (Figure 1) 1]. This process may 
create sufficiently high pressures in the CO cores to nu- 
cleate diamond crystallites @ which continue to grow 
under irradiation [6j. Studies of irradiated CO with en- 
capsulated metal crystallites provide additional evidence 
for large pressures which lead to a substantial shift in 
the melting temperature of the encapsulated material 0] ■ 
Hydrostatic pressures above 20GPa were demonstrated 
by determining the lattice constant of Au nanoparticles 
encapsulated in irradiated CO Q. Pressurization also oc- 
curs in MWNT. frradiation-induced MWNT contraction 
was observed to reduce the lattice constants of encap- 
sulated metals and carbides, indicating pressures up to 
40 GPa, and to induce plastic deformation at stresses 
close to the theoretical shear stress 0]. Shell shrink- 
age has been attributed to bond re-arrangement around 
irradiation-induced point defects. Simulations indicate 
reconstruction of di-vacancies which leads to saturation 
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FIG. 1: CO during 1.2MeV e" irradiation at 1000 K; coher- 
ent shells with decreased spacing towards the core 



of dangling bonds and concomitant 'loss' of unoccupied 
sites in the sp 2 network [HJ [H| ■ In single- walled nan- 
otubes this decreases the average nanotube radius by an 
amount that is proportional to the di- vacancy concentra- 
tion [9(. In multi-shell structures such as CO and MWNT 
the situation is more complex owing to the migration and 
reactions of point defects within multiple shells [!, [f| ■ 
Existing rate equation models for the point-defect dy- 
namics in irradiated CO [13[ do not account for the loss of 
lattice sites by di- vacancy re- arrangement. In the present 
work, we formulate a model of the dynamics of point- 
defects in closed-shell graphene nanostructures which for 
the first time accounts for the di-vacancy reconstruction 
mechanism, the concomitant geometrical changes, and 
the evolution of internal stresses. 

We consider a spherical or cylindrical graphene nanos- 
tructure consisting of M concentric shells Si, 1 < i < M. 
The outermost shell is labeled Si . Ni is the number of sp 2 
bonded sites on Si. (Note that the shells may not pos- 
sess perfect Graphene structure but contain pentagons, 
Stone- Wales defects [l4j], or more complex bonding pat- 
terns.) The stress-free radius of Si is r t — [A r i yl*/(47r)] 1 / 2 
for CO and n = NiA*/(nl) for MWNT where A* is the 
average area per atom which we approximate by the value 
for planar graphene and I 3> the MWNT length. The 
spacing between shells i and i + 1 is di := rj — ri+j.. 
The number of vacant lattice sites on Si is nj , and the 
number of interstitials contained between Si and 5^+1 is 
n\. These interstitials diffuse two-dimensionally on an 
'interstitial shell' of area A; (A; = 7r(r.; + r^+i) 2 for CO). 
Diffusion of vacancies is considered negligible because of 
their significantly higher migration energy [l2j]. 

We envisage the following processes: (i) Irradiation 
displaces atoms at rate <!> from sp 2 bonded sites. Dis- 
placement of an atom from Si creates a vacancy on this 
shell and, with equal probability, an interstitial either be- 
tween Si and iSi+i {n\ — > n\ + 1) or between Si and Si-i 
(n?_j — > n\_ 1 + 1). Accounting for vacant lattice sites, 
the total rate for this process is &(Ni — nj). (ii) Recom- 
bination of interstitials with vacancies on adjacent shells 
is modelled as a bimolecular diffusion-controlled reaction 
where the reaction constant K « A*v$ exp[— E^/ (k-sT)] 
approximately equals the diffusion rate of interstitials be- 
tween their enclosing shells, vq « 10 13 s _1 is a frequency 
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of the order of the Debye frequency and Em the inter- 
stitial migration enthalpy. Stress gradients affect the re- 
combination rates: If there is exactly one vacancy on Si 
and one on <Si+i, then in absence of internal stresses an 
interstitial diffusing between these shells has equal prob- 
abilities Pi.i — Pi,i+i = 1/2 to recombine with cither 
vacancy. If the graphene shells experience different in- 
plane stresses Cj, on the other hand, we assume that this 
leads to preferential recombination of interstitials in the 
direction of increased tensile/reduced compressive stress. 
The modified recombination probabilities follow from 
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where V — 8E W /da characterizes the dependence of 
the vacancy energy E in graphene on in-plane tensile 
or compressive stresses. Thus, creation and recombina- 
tion of interstitials may lead to a net flux of matter in 
the direction of decreasing compressive/increasing ten- 
sile stresses, (iii) Shells shrink when irradiation creates 
di-vacancies which reconstruct to saturate the dangling 
bonds, resulting in two pentagons and one octagon in 
the bond network (see e.g. [l0|, HH). The resulting de- 
fects do not recombine with interstitials as the energy of 
a reconstructed divacancy is less than that of a single va- 
cancy 15[ - we may thus think of di-vacancy formation 
on Si as the removal of two sites from the sp 2 network, 
Ni — > Ni — 2, and of two vacancies, nj — > nj — 2. For sim- 
plicity, we assume that A* is not changed by this process 
which thus reduces the shell surface by 2 A* and the mean 
shell radius by A* /(4wri) for CO. With vacancy diffu- 
sion assumed negligible, di-vacancy creation is governed 
by the displacement of nearest neighbors of pre-existing 
vacancies, and occurs at rate 3<&nJ . 

Combining processes (i) to (iii), we arrive at coupled 
rate equations for the vacancy and interstitial numbers: 



d t nj 



$(Ni - 7nJ) 



Kn) 



n 



i 

(-1 



d t n\ = <S>(N t + N i+i - nj) - nY +1 )/2 



Pi-i 
V\ „v 



Kn 



P H i _L p '+ 1 
1 1,1 . T -i»,t+l . 



(2) 



(3) 



These equations need to be modified in the core (i = M ) 
where Pm.m = 1 an d N^+i = n Ji+i = 0, and on the 
outermost shell S\ where Po.i = as atoms that are 
displaced outwards do not recombine with vacancies but 
leave the system. The rate of contraction of Si due to 
di-vacancy creation is given by dtri = — 3^nf A* / (Awri) 
for CO. 

To obtain a closed description, the geometry needs to 
be related to the internal stresses acting in the struc- 
ture. To this end, we consider the inner and outer pres- 
sure acting on Si because of Van der Waals interactions 
with the neighbouring shells. If there is a difference 



between inner and outer pressure, force equilibrium re- 
quires this to be balanced by a hoop stress. Consider- 
ing Si as a thin shell of radius r, and thickness di, the 
hoop stress is given by = (/,: — / (adi) where 

fi is the van der Waals force per unit area acting on 
the shell from inside due to its interaction with shell 
i + 1, a = 1 for MWNT and a = 2 for CO. To eval- 
uate the forces, we refer to experimental data for pla- 
nar graphite. In a third-order polynomial approxima- 
tion, the force per unit area acting on Si from inside 
is approximated by /; = c 33 [e£ + K 1 (e c i ) 2 + K 2 [{e^) 3 } 
where el — (di/do) — 1 and we use for the elastic con- 
stant C33 the value for bulk graphite, C33 ~ 36 GPa [l6j]. 
do ~ 0.335nm is the shell spacing in absence of internal 
stresses which we take to equal the c-axis lattice spacing 
of bulk graphite. The polynomial coefficients are deter- 
mined by fitting to experimental pressure-compression 
data [nf, yielding K x w 0, K 2 « 130, i.e., the third- 
order non-linear corrections prevail; they are substantial 
even for moderate degrees of compression. The bound- 
ary condition at the outer surface of the nanostructure is 
po,i = 0. To evaluate di and Ti, we consider the shells as 
rigid and approximate the shell radii by their stress-free 
values. This is justified by the extreme elastic anisotropy 
of graphitic structures, as the in-plane elastic modulus of 
the graphene layers exceed the modulus in the perpen- 
dicular direction by a factor of 30. 

With an interstitial migration energy Ef 1 — 0.8 eV 
[l2l ] and typical irradiation rates of $ « 10~ 4 ... 1 dpa/s, 
the model is applicable to temperatures T > 230 K— 310 
K, below which defect concentrations become so large 
that the shells lose coherency. We thus exclusively con- 
sider irradiation at above-ambient temperatures and ex- 
plore the experimentally accessible range of temperatures 
and irradiation rates, 400 K < T < 1200 K and 10~ 4 
s _1 < $ < 1 s _1 . The relaxation volume V which gov- 
erns the stress dependence of the vacancy energy may be 
estimated from the structural relaxation of atoms around 
a vacancy, indicating V ~ 0.14A*do [18]. As initial con- 
dition we consider a stress- and defect free structure, i.e., 
n\ = nj = and di — do for all i. The CO radius is 
R(t) = r\{t) and the initial radius is Rq — Ndo. 

The evolution of an irradiated CO proceeds in two 
stages, characterized by stress build-up followed by dia- 
mond nucleation and growth or, if pressures are insuffi- 
cient for diamond nucleation, by self-similar contraction 
of the nanostructure. At the onset of irradiation, quasi- 
stationary point defect populations build up. Above 
room temperature these are low everywhere with the ex- 
ception of Si where loss of atoms by sputtering leads to a 
gradually increasing excess of vacancies. Displacements 
of atoms next to these vacancies lead through di-vacancy 
arrangements to gradual shrinkage of Si which then ex- 
erts pressure on 1S2 while itself experiences a tensile hoop 
stress. These stresses modify the point defect recombi- 
nation rates: Interstitials created between the two shells 
now preferentially recombine with vacancies on Si. Ac- 
cordingly, the shrinkage rate of Si is reduced while that 
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FIG. 2: Evolution of the vacancy number n( (multiplied with do/n), hoop stress Oh, recombination probability Pi_H, and shell 
spacing di during the first stage of irradiation of a CO with 16 shells; <E> = 0.01 dpa/s, T = 570if, other parameters see text. 
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FIG. 3: Stage of constant-rate contraction of an irradi- 
ated CO, parameters as in Figure [2] 



of 1S2 is increased as excess vacancies build up. This re- 
lieves the compressive stress on 52 and transfers it to 1S3 
while the spacing between S\ and 1S2 henceforth remains 
approximately constant. The process then repeats itself 
with 54, and so on. As a result, during the stage of stress 
build-up the onion consists of a set of outer shells that 
shrink at approximately equal rate. These outer shells 
are under tensile stress and have high vacancy concentra- 
tions. They surround a stress-free non-contracting core 



with low vacancy concentration. The outer shells are 
separated from the core by an interface where large com- 
pressive stresses reduce the vacancy recombination rate 
and thus build-up excess vacancies and induce shrinkage. 
The evolution of the nanostructure is characterized by 
inwards propagation of the interface and the accompa- 
nying pressure wave. The process is illustrated in Figure 

m 

The height of the compressive stress peak at the inter- 
face between the contracting outer layers and the station- 
ary core increases as the interface propagates inwards. 
Systematic simulations show that the compressive stress 
maximum (the highest stress acting on an inner layer) 
reached during onion contraction is approximately 

~ r NkBT U\ 

^c,max ~ J ^ \?*) 

where / 1.7. For an onion with N = 10 deforming 
at T = 1000K, compressive stresses may reach up to 27 
GPa. The further evolution of the nanostructure depends 
on whether the peak compressive stress is sufficient for 
diamond nucleation. If this is not the case, we reach 
a stage where all shells contract at approximately equal 
rate proportional to $ (Figure [3]) and CO contraction 
goes along with the progressive 'disappearance' of shells 
in the onion centre. The shell spacing profile is charac- 
terized by a single scaling function d[r / R) which depends 
only on the ratio k^T/V and remains invariant as R de- 
creases - we thus speak of self-similar contraction. 
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FIG. 4: Top: shell spacing profiles during self-similar 
contraction of CO with different initial shell numbers N 
and irradiation temperatures T. Empty symbols: R/Ro 
= 0.9, cross center symbols: R/Ro = 0.6, full symbols: 
R/Ro = 0.3. Bottom: Contraction of a CO with N = 14 
during irradiation at 1000 K with 1.2 MeV e~ (flux den- 
sity 150 A/cm 2 ) [13| ; the experimental data were taken 
after irradiation doses of 6dpa, 12dpa and 60dpa; Cross- 
center symbols: CO with 28 shells after irradiation under 
the same conditions to approximately 30 dpa. 



We compare the predictions for this contraction 
regime with observations of irradiated CO by Banhart 
and Lyutovic USE]- I n these studies, irradiation was 
carried out in a 1.2 MeV high- voltage electron microscope 
with typical beam intensities of 4.5 x 10 5 A/m 2 , corre- 
sponding to typical displacement rates of 10~2 dpa/s. 
Irradiation at temperatures around 1000K induced pro- 
nounced compression after doses of a few dpa [19j . Shell 
radii were determined from electron micrographs by trac- 
ing the contrast peak which indicates the location of a 
shell and determining the enclosed area. Shell spacing 
profiles obtained at different irradiation doses are shown 
in Figure |4] for onions with initially 14 and 28 shells (the 
28-shell onion is shown in Figure 1). The figure demon- 
strates that the onions indeed contract in a quasistation- 
ary manner with approximately time-independent shell 
spacings. The calculated shell spacing profiles are in good 
agreement with the experimental data. 

In CO above a critical size and/or irradiation temper- 
ature, the calculated pressures are sufficient to explain 
the nucleation of diamond in the CO core. Computations 
indicate a critical stress of about 16 GPa for sp 3 bond 
formation during compression of MWNT and 12 GPa 
for Graphite [2(|, which allows us to estimate the 'phase 
boundary' for diamond nucleation in CO as NT[K] I 
4500-6000. Once nucleated, the continuous conversion of 
sp 2 bonded graphene to sp 3 bonded diamond maintains 
the pressure at the shell-core interface at the critical 
level required for diamond growth. A study of the 
growth regime needs, however, to take into account 
the simultaneous operation of other irradiation-driven 
diamond growth mechanisms @ and is thus beyond the 
scope of the present letter. 



[1] F. Banhart and P. M. Ajayan, Nature 382, 433 (1996). 

[2] F. Banhart, J. Appl. Phys. 81, 3440 (1996). 

[3] F. Banhart, T. Fuller, Ph. Redlich, and P. M. Ajayan, 

Chem. Phys. Letters 269, 349, (1998). 
[4] P. Wesolowski et. al., Appl. Phys. Letters 71, 1948 

(1997). 

[5] Ph. Redlich, F. Banhart, Y. Lyutovich and P. M. Ajayan, 

Carbon 36, 561 (1998). 
[6] M. Zaiser and F. Banhart, Phys. Rev. Letters 79, 3680 

(1997). 

[7] F. Banhart, E. Hernandez, and M. Terrones, Phys. Rev. 

Letters 90, 185502 (2003) 
[8] L. Sun et. al., Phys. Rev. Letters 101, 156101 (2008) 
[9] L. Sun et. al, Science 312, 1199 (2006). 
[10] K. R. Bates and R. E. Scuseria, Theor. Chem. Acc. 99, 

29 (1998). 

[11] M. Sammalkorpi et. al., Phys. Rev. B 70, 245416 (2004). 



[12] F. Banhart, Rep. Progr. Phys. 62, 1181 (1999). 
[13] W. Sigle and Ph. Redlich, Phil. Mag. Letters 76, 125 
(1997). 

[14] A. J. Stone and D. J. Wales, Chem. Phys. Letters 128, 
501 (1986). 

[15] A. V. Krasheninnikov, Chem. Phys. Letters 418, 132 
(2006). 

[16] K.H. Michel and B. Verberck, phys. stat. sol. (b) 245 

21772180 (2008). 
[17] R. W. Lynch and H. G. Drickamer, J. Chem. Phys. 44, 

181 (1966). 

[18] A. P. P. Nicholson and D. J. Bacon, phys. stat. sol. 28 
613 (1975). 

[19] Y. Lyutovich, thesis, Stuttgart University 2000. 

[20] W. L. Guo et. al., Phys. Rev. Letters 93, 245502 (2004). 



